arXiv:1501.05638vl [cond-mat.str-el] 22 Jan 2015 


A Measure of Monopole Inertia in the Quantum Spin Ice Yb 2 Ti 207 


LiDong Pan 1 , N. J. Laurita 1 , Kate A. Ross 1 ’ 2 , Edwin Kermarrec 3 , Bruce D. Gaulin 3 ’ 4 ’ 5 , and N. P. Armitage 1 

1 Institute for Quantum Matter, Department of Physics and Astronomy, Johns Hopkins University, Baltimore, 
Maryland 21218, USA 

2 NIST Center for Neutron Research, NIST, Gaithersburg, Maryland 20899, USA 

3 Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 4M1, Canada 

4 Brockhouse Institute for Materials Research, McMaster University, Hamilton, Ontario, L8S 4M1, Canada 

5 Canadian Institute for Advanced Research, 180 Dundas St. W, Toronto, ON, M5G 1Z8, Canada 

(Dated: January 23, 2015) 


An important and continuing theme of modern 
solid state physics is the realization of exotic exci¬ 
tations in materials (e.g. quasiparticles) that have 
no analogy (or have not yet been observed) in the 
actual physical vacuum of free space. Although 
they are not fundamental particles, such quasi¬ 
particles do constitute the most basic description 
of the excited states of the ”vacuum” in which 
they reside dU [2} [3]) . In this regard the magnetic 
textures of the excited states of spin ices , mag¬ 
netic pyrochlore oxides with dominant Ising inter¬ 
actions, are proposed to be modeled as effective 
magnetic charge monopoles (j4j[5]). Recent inelas¬ 
tic neutron scattering experiments have estab¬ 
lished the pyrochlore material Yb 2 Ti 2 07 (YbTO) 
as a quantum spin ice, where in addition to the 
Ising interactions there are substantial transverse 
terms that may induce quantum dynamics and - 
in principle - coherent monopole motion (|6j [7ft . 
Here we report a combined time domain terahertz 
spectroscopy (TDTS) and microwave cavity study 
of YbTO to probe its complex dynamic magnetic 
susceptibility. We find that the form of the sus¬ 
ceptibility is consistent with monopole motion 
and a magnetic monopole conductivity can be de¬ 
fined and measured. Using the unique phase sen¬ 
sitive capabilities of these techniques, we observe 
a sign change in the reactive part of the magnetic 
response. In generic models of monopole motion 
this is only possible through introducing inertial 
effects, e.g. a mass dependent term, to the equa¬ 
tions of motion. Analogous to conventional elec¬ 
tric charge systems, measurement of the conduc¬ 
tivity’s spectral weight allows us to derive a value 
for the magnetic monopole mass, which we find 
to be approximately 1800 electron masses. Our 
results establish the magnetic monopoles of quan¬ 
tum spin ice as true coherently propagating quasi¬ 
particles of this system. 

Quantum spin ice has received considerable recent at¬ 
tention in the search for quantum spin liquids, as a pos¬ 
sible realization for this novel state of matter in which 
highly quantum entangled spin degrees of freedom evade 


conventional long range order down to the lowest temper¬ 
atures ©■ Magnetic pyrochlore oxides, in which mag¬ 
netic ions sit at the vertices of corner sharing tetrahedra, 
provide a three dimensional system where spin ice states 
are found, when provided with appropriate spin interac¬ 
tions and anisotropy ©■ In classical spin ice materials, 
such as Dy 2 Ti 2 07 and Ho 2 Ti 2 07 , large classical spins 
are forced by strong crystal field anisotropy in the local 
<111> direction, with a primarily ferromagnetic Ising 
spin interaction. The resulting ground state obeys the 
Bernal-Fowler ice rules, where each tetrahedron adopts 
the so-called “two-in, two-out” configuration. This is 
equivalent to the proton configurations in water ice (two- 
close, two-far) and hence the classical spin ices are found 
to similarly possess an extensive low temperature residual 
entropy. Flipping a single spin (e.g, a dipole excitation) 
in the spin ice creates a pair of magnetic monopoles in 
neighboring tetrahedra, which can then subsequently be 
separated by additional spin flips resulting in deconfined 
monopoles OED. 

The slow dynamics of the monopoles in the classical 
spin ices are still a subject of investigation, but they are 
believed to be driven by the strong fluctuating transverse 
component of the dipolar field arising from thermally 
fluctuating neighboring spins©. In contrast, in mate¬ 
rials like YbTO, with the addition of finite transverse 
terms in the spin Hamiltonian, monopole dynamics be¬ 
come inherently quantum and change the situation dra¬ 
matically (Ushiihiid. The exchange interaction param¬ 
eters for the magnetic pyrochlore oxide YbTO have been 
obtained by recent inelastic neutron scattering (INS) at 
high magnetic fields and later confirmed by other 
experimental and numerical studies 0 03} H3D • While 
the largest interaction is Ising and ferromagnetic, sig¬ 
nificant non-Ising contributions exist. The crystal field 
structure of the material allows the well isolated ground 
state Kramers doublet to be treated as an effective spin 
1/2 moment in the low energy sector (IT51 fTO[) . These re¬ 
sults establish YbTO as an exchange quantum spin ice. 
Although the exact nature of its ground state is under 
debate © [13} 17; T5} [12} HU; 211 122), at temperatures 
above a first order transition (T c ~ 260 mK), YbTO is 
believed to be in a spin ice-like state, but one whose dy¬ 
namics are determined by non-Ising terms, which allow 
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FIG. 1 Schematic of spin ice and experimental setup Schematic illustration of a quantum spin ice with deconfined monopole 
and antimonopoles (shown as red and blue spheres), together with incident and transmitted terahertz pulses. In the pyrochlore 
lattice, magnetic rare earth ions sit at the vertices of corner-sharing tetrahedra, which are shown as grey spheres in the figure. 


the quantum tunneling of magnetic monopoles between 
tetrahedra. Despite the advances in understanding its 
behavior, there has been as of yet no definitive evidence 
for spin ice like correlations in YbTO. With quantum 
spin liquids lacking an experimentally verified smoking 
gun, the characterization of the magnetic monopoles in 
the quantum spin ice state is a relevant and urgent task. 

Here we study the dynamics of the magnetic response 
in the quantum spin ice regime of YbTO with combined 
TDTS (schematic illustrated in Fig. 1) and microwave 
cavity techniques. Experiments were performed down to 
a temperature of 1.4 K, which is well below the peak in 
the low temperature magnetic heat capacity and within 
the purported quantum spin ice regime (UJ) (and well be¬ 
low the estimated scale of the mean-field T c ~ 3.4 K 
©)■ Low frequency electromagnetic spectroscopies can 
provide highly detailed and accurate information of the 
dynamic response of the material in the relevant energy 
range. With the recent advances of phase sensitive op¬ 
tical techniques such as the ones employed here, com¬ 
plex material response function can be directly obtained 
which brings crucial knowledge to our analysis. The mag¬ 
netic response of this material is consistent with a model 
of massive monopole motion. 

The physical quantity we directly measure in a TDTS 
experiment is the complex transmission function t(ui). 
Frequency and temperature dependence of the magni¬ 
tude of t(u>) measured from TDTS in zero field are shown 
in the intensity plot in Fig. 2(a). The spectra are nor¬ 
malized by a scan at 50 K, a temperature well above the 
onset of magnetic correlations. An increase in the low 
energy dissipation is observed as temperature decreases 
which is unusual for a large gap insulator. Such an ef¬ 
fect is unlikely to come from lattice effects, as the lowest 
infrared active optical phonon in YbTO has an energy 
of 2.25 THz. Moreover, there are no reported dielectric 
anomalies at low temperature US USD- The feature is 
most pronounced for frequencies below 0.4 THz, and for 


temperatures below 10 K. This trend is reminiscent of 
the low frequency magnetic susceptibility, as this temper¬ 
ature range coincides with the onset of spin correlations 
in YbTO and the eventual crossover to the quantum spin 
ice state. The energy scale of this feature matches the 
scale of the spin interactions obtained from INS ©. 

To further demonstrate the magnetic nature of the 
low frequency absorption, we show the magneto-optical 
TDTS spectra of YbTO in the Voigt geometry in Fig. 
2b. The data shows the THz transmission amplitude as 
a function of frequency and applied magnetic field taken 
at 1.6 K, normalized by the spectrum at 7 T. The field 
dependence of the absorption features proves the mag¬ 
netic origin of the zero field dissipative mode. The be¬ 
havior of this system in magnetic field have been studied 
in detail recently (jH [T4j HZ|). At high fields, YbTO is 
in a spin polarized state, where excitations are magnon 
and two-magnon excitations (two dark lines diagonally 
crossing Fig. 2(b)). As the field decreases, those dipole 
excitations crossover to weakly confined quantum string¬ 
like excitations that connect monopole - antimonopole 
pairs, which is signified by their increasing effective g- 
factors below 3 T (O. Eventually the magnetic reso¬ 
nance modes lose their sharp structure and the spectra 
exhibit a broad, diffuse feature at zero field which is con¬ 
sistent with the magnetic monopoles being deconfined in 
the zero field quantum spin ice state. We show below 
that, despite its apparently broad profile, this absorp¬ 
tion exhibits features that can be connected to coherent 
monopole transport. 

Given the separation of frequency and temperature 
scales of the magnetic and lattice degrees of freedom, 
we can extract the complex magnetic response from our 
complex transmission. The spectra at low temperatures 
are referenced to a scan taken above the onset tempera¬ 
ture of the magnetic correlations. More specifically, we 
have In [t(T)/t(T re f)\ = in'uJx(oj)d/2c, here x(w) is the 
complex dynamic magnetic susceptibility of the system 
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FIG. 2 Transmission as a function of temperature and field (a) The intensity plot of the transmission amplitude from TDTS 
measured in zero field, as a function of frequency and temperature, normalized with the spectrum at 50 K. (b) Intensity plot 
of the transmission amplitude as a function of frequency and field, normalized with the spectrum at 7 T. The data in this plot 
are measured at 1.6 K with H oc || Hd c T k thz (Voigt geometry). 


in the limit q —>■ 0, c is the speed of light, and d is the 
sample thickness. We make the reasonable assumption 
for an insulator at low temperature and frequency that 
the index of refraction n' from dielectric effects is real, 
with little temperature and frequency dependence below 
20 K in the THz range. Further details of the analysis 
are included in the supplementary information (SI). 

In analogy with conventional electric charge conductiv¬ 
ity, we can define a magnetic monopole conductivity as 
k = ujx(u))/i. In Fig. 3(a) and (b) we show the real and 
imaginary parts of n as a function of frequency at sev¬ 
eral temperatures. Data from TDTS are shown as open 
symbols. Other than the increasing dissipative response 
at low temperature, the most striking feature of the data 
is the trend of the frequency dependence of the reactive 
response (k"), suggesting a sign change of x' a t a fre¬ 
quency just below the low frequency cut off for TDTS. 
This sign change is further confirmed by the addition of 
microwave cavity measurements at 18.5 GHz, which are 
shown in Fig. 3 as solid dots. The combined zero field 
measurements point to the existence of a zero crossing of 
X’ at a frequency slightly lower than 0.1 THz. In the inset 
to Fig. 3, we show the corresponding complex magnetic 
susceptibility at low temperatures. Also, an external dc 
magnetic field pushes the zero crossing to higher frequen¬ 
cies where reliable TDTS data can be achieved, provid¬ 
ing further support for this feature. Details of the cavity 
data analysis and these magneto-optical TDTS results 
are included in the SI. 

The existence of the sign change in the reactive mag¬ 
netic response function is a remarkable feature by itself. 
To describe the dynamics of classical spin ices, Ryzhkin 
proposed a Debye-like model with an equation of motion 
written as J; = pn ljn (q l)rn Tl — ^■FJl) (j4j). Here J is 
the monopole current, p is the monopole mobility, n m is 
the monopole density, and q m is the monopole charge. 
rji is an index that takes on values ±1 for the charge 


polarity with the subscript i denoting the monopole po¬ 
larity. The second term in the Ryzhkin model describes 
a reaction field which originates as a consequence of the 
configuration entropy due to the motion of monopoles, 
which prohibits a true direct monopole current even in 
the thermodynamic limit. Here i”2 = M/g m is a configu¬ 
ration vector that is proportional to the system’s magne¬ 
tization while $ = -^aksT is a temperature dependent 
proportionality constant. The Ryzhkin model describes 
essentially overdamped monopole motion, and generates 
a frequency dependent susceptibility that is Debye-like in 
form. The model has enjoyed great success in the study 
of classical spin ice materials (El, but fails completely in 
the attempt to account for the sign change observed in k" 
in the quantum spin ice YbTO. Moreover, the Ryzhkin 
model predicts flat frequency dependence of k' whereas 
the data decreases at high ui. The only way to capture 
the sign change in k" and the fall off at high uj in k! 
in any classical model of monopole transport is by the 
inclusion of the inertial effects i.e. a term which is pro¬ 
portional to the second time derivative of the coordinate 
of interest. This can be seen from the fact that the lead¬ 
ing dependence of the susceptibility in powers of l/(— iui) 
is demonstrative of the dominant term in the equations 
of motion at a particular frequency. 

The inclusion of a phenomenological term describing 
inertial effects in the equation of motion for magnetic 
monopoles may extend the Ryzhkin model to the quan¬ 
tum spin ice systems. Our equation of motion then 
reads J, = pni^ m {qi ^ m H — ^$17) — Jj/ 7 , where 7 is a 
monopole current relaxation rate. From standard defini¬ 
tions, p = qm/ r T m m- Here m m is an effective monopole 
inertial mass that arises in the low energy sector through 
non-Ising exchange that lead to monopole tunneling. 
There is no more contradiction inherent in using a clas¬ 
sical model to describe a quantum spin ice than there 
is to use a classical model to describe the conduction of 
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FIG. 3 Real and imaginary parts of measured k with fitting Real (a) and imaginary (b) part of the magnetic conductivity 
k = a }\(uj)/i as a function of frequency at several temperatures. Solid symbols are data obtained with microwave cavity 
measurements, while open symbols show data from TDTS. Solid lines are fitting curves as described in the main text. Inset 
shows the corresponding data and fitting of the complex magnetic susceptibility at two temperatures. 


thermally excited charge in a semiconductor where clas¬ 
sical inertia also arises through inherently quantum tun¬ 
neling. Further motivation for this approach is given in 
the SI. Solving in a similar fashion as Ryzhkin yields a 
Drucle-Lorentz style expression for the complex dynami¬ 
cal susceptibility x( w ) 


xM 


1 — iojT — w 2 t /7 


(1) 


The inertial effects of magnetic monopoles may arise 
as a consequence of their ability to coherently tunnel be¬ 
tween tetrahedra through the finite non-Ising terms in 
the spin Hamiltonian. The quantum spin ice system is 
distinguished by its much lower current damping rate 
than the classical spin ice systems. k(oj) can be fitted 
to Eq. [TJ which contains three independent parame¬ 
ters, xo = 7) and 1/r = n m /r$/g m . Here xo 

is the dc value of the magnetic susceptibility, while 1 /r 
is the magnetization relaxation rate that also appears 
in Ryzkhin’s theory. 7 is a new rate introduced by our 
treatment that parametrizes monopole momentum relax¬ 
ation. Fits shown in Fig. 3 as solid lines are obtained by 
simultaneously fitting the real and imaginary parts of 
k(ui) including both the microwave cavity and the TDTS 
data. Although the fits are not perfect, it is important to 
emphasize how restrictive the fitting procedure is (even 
despite the frequency gap between the TDTS and mi¬ 
crowave data) due to the Kramers-Kronig relation be¬ 
tween the complex values of k(w). k 1 at low frequency 
cannot be changed without affecting k" at higher fre¬ 
quency and vice versa. This simultaneous fitting of the 
real and imaginary parts of the data ensures that the 
essential physics is well captured and constrained. The 
ability of this model to capture the salient features pro¬ 
vides strong support for the phenomenological descrip¬ 
tion of quantum spin ice developed here, although a more 
sophisticated model that incorporates a distribution of 
relaxation times as done recently in the classical spin 


ices (pH) might potentially improve the fits. As the spec¬ 
tra themselves evolve continuously with temperature, we 
have performed these fits up to a reasonably high tem¬ 
perature of 10 K, which is presumably well out of the 
regime of any spin ice correlations. As discussed be¬ 
low, in the paramagnetic regime at temperature above 
the mean-field T c , we believe that fits to Eq. 1 should 
be considered only phenomenological without any corre¬ 
spondence of the fit parameters to monopole character¬ 
istics. 

The temperature dependence of the fitting parameters 
is shown in Fig. 4. The excellent agreement between 
Xo obtained from the current experiment and the val¬ 
ues reported in Ref [T3] with magnetization measurement 
lends further support to our phenomenological model 
of the monopole dynamics in quantum spin ice at low 
temperatures. The satisfying match between the xo ob¬ 
tained from our measurements and the dc value also sug¬ 
gest (due to the above mentioned Kramers-Kronig con¬ 
straints) there are no significant features at even lower 
frequency in the magnetic susceptibility of YbTO. Both 
relaxation rates 7 and 1 /r first show a slow decreasing 
trend upon cooling down to ~ 4 K (of order the mean- 
field T c ©), with a faster dependence below this temper¬ 
ature, before apparent saturation below 2K and down to 
our lowest measured temperatures. We believe the data 
can be interpreted in terms of monopole motion below the 
mean-field temperature scale. The magnetization relax¬ 
ation rate we find at low temperature (1/r « 230 GHz) is 
consistent with the limits put on it by neutron spin echo 
experiments {l/T ec ho > 250 GHz) (f?7| 1 . At this time, 
the source of the temperature dependence of the momen¬ 
tum relaxation 7 is unclear. But in analogy with charge 
systems, its possible that the temperature independent 
offset (e.g. a residual resistivity) is caused in part by im¬ 
purities. The present system has defects at the 1% level 
that are describable in terms of ’’weak stuffing” ED, and 
possible oxygen non-stoichiometry. 

The inclusion of inertial effects also repairs a defect 
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Temperature (K) 

FIG. 4 Temperature dependence of the fitting parameters 
for the extended Ryzhkin model describe in the main text. 

(a) xo obtained from the current experiment (dots), and from 
dc susceptibility measurements reported in Ref. (|13[) (line). 

(b) Relaxation rates 7 , r _1 . In this plot, rates are divided 
by 2n to put in units to compare to spectral plots. To com¬ 
pare these quantities to relaxation in time, quantities should 
have the 2n removed, (c) Magnetic spectral weight (plasma 
frequency squared) obtained from fitting. 


of Ryzhkin’s Debye-like model in that it is not, strictly 
speaking, a fully mathematically consistent magnetic re¬ 
sponse function as it does not fall off fast enough at 
high 10 to satisfy the zero temperature first moment sum 
rule / 0 °° duj wy" Q (q,w) = gp([Sq, [^, S-q]]) for lo¬ 
cal spins (1251 I22|) . Here S'" = ]T\ e - jq R j S'" j s the 
Fourier transform of lattice spin operators and H. is the 
spin hamiltonian. A similar violation of a sum rule is 
a well-known defect of the Debye model in its applica¬ 
tion to electronic charge relaxation in dielectrics ®. 
In general the first moment of the electric susceptibility 
must obey the relation dui = /o°° ^ a ' = f 
where w 2 = n e q 2 /m e defines the plasma frequency ui p , 
a is the electric charge conductivity, n e is the electron 


density and m e the electron mass. Debye-like models do 
not fall off fast enough at high frequency to satisfy this 
sum rule. The addition of inertial terms to the equations 
of motion in Debye-like models gives \ a high frequency 
asymptote that goes like 1 /(—iw) 2 which insures the first 
moment’s integrability flU}. Although at low frequencies 
in strongly dissipative media inertia can be neglected, at 
high enough frequencies it must become relevant to sat¬ 
isfy the sum rule. In electronic charge systems frequent 
use is made of this sum rule to determine the mass of 
charge carriers (if the charge density is known). 

With the mapping of the form of the monopole con¬ 
ductivity described here to the Drude-Lorentz model, we 
expect a similar sum rule to exist for the magnetic spec¬ 
tral weight. To complete the analogy in the present case, 
we can extract a spectral weight from the fits that is 
related to a monopole plasma frequency and more fun¬ 
damental monopole parameters as w 2 = n m qm/m m . The 
temperature dependence of the spectral weight is shown 
in Fig. 4(c), which illustrates the decreasing trend of the 
magnetic spectral weight upon warming. Further analy¬ 
sis of the fits to the susceptibility of the extended Ryzhkin 
model allows us to extract an effective mass in m of the 
monopoles from the fitting parameters. Here we eval¬ 
uate the low temperature value of m m = n m q^ n T/xo r Y- 
We use the monopole density n m from the results of 
Ref. [7J which used the Numerical Linked Cluster method 
to successfully model the heat capacity of YbTO. Us¬ 
ing this computed value for n mi we determine a value of 
~ 1.6 x 10~ 27 kg or ^1800 m e for the monopole mass at 
temperatures below 4 K (for further details see SI). 

As discussed earlier, the mass of the monopoles orig¬ 
inates from their ability to tunnel between sites in the 
quantum spin ice systems as a consequence of the finite 
non-Ising exchange terms. Unfortunately, reliable theo¬ 
retical estimates for the inertial mass of the monopoles 
in YbTO in this temperature range do not exist. How¬ 
ever, it is believed that their tunneling rate is governed 
by the transverse exchange term J z ± 0®. We obtain 
a rough estimate for the monopole mass with a tight- 
binding model of charges on the diamond lattice (the dual 
lattice of spins). In the q —> 0 limit a calculation gives 
to ~ 4fi 2 /(a J z± d 2 ) where d is the diamond lattice unit 
cell parameter and a is a temperature dependent dimen¬ 
sionless constant of order unity. With a = 1, such a treat¬ 
ment gives an effective monopole mass of ~ 2.0 x 10 -27 
kg (~ 2200 TO e ), which corresponds closely to our experi¬ 
mental value. As noted above, our spectra evolve contin¬ 
uously into the high temperature paramagnetic regime. 
Although one can continue to fit the data to the form 
of Eq. 1, the fitting parameters become increasingly un¬ 
physical above the temperature scale of the mean-field 
T c . At 10K one would determine the fitted mass is 10,000 
m e and strongly increasing (see SI), which we do not be¬ 
lieve has any physical significance. This agreement of the 
fitting parameters with a monopole model at low tem¬ 
peratures and the disagreement at higher temperatures 
is further testament to the validity of our interpretation. 
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Following the mass determination of the monopoles 
and mapping to monopole conductivity, a mobility can 
also be obtained from an analog of the standard expres¬ 
sion /i = qm/'inj, which yields a weakly temperature 
dependent mobility of ~ 100 ms -1 T -1 for the lowest 
temperature of our measurement at 1.4 K (a detailed 
plot can be found in the SI). This value is several orders 
of magnitude larger than the mobility found in classical 
spin ice materials (PH). This, along with the tempera¬ 
ture independent regime of the relaxation rates at low 
temperature, points to a clear distinction between the 
monopole transport in the two classes of spin ice mate¬ 
rials. This supports the idea that although in principle 
the mass would get larger (and inertial effects naively 
would increase) in the classical Ising limit, quantum tun¬ 
neling from transverse exchange terms is eventually su¬ 
perseded by incoherent monopole hopping. We note that 
our TDTS experiments of Dy 2 Ti 2 C >7 show no temper¬ 
ature or field dependent features in the spectral range 
reported here. The demonstration of the sign change in 
the reactive part of the complex magnetic response func¬ 
tion is strong evidence of the monopoles’ inertial effects 
in quantum spin ice. With this, their identity as the true 
coherent quasiparticles of YbTO is further established. 


I. METHODS 

TDTS experiments are performed in a home built spec¬ 
trometer with two dipole switches as THz emitter and 
detector. For experiments in magnetic field, a cryogen 
free superconducting magnet with optical access is used 
for hosting the sample. Single crystal YbTO samples are 
grown with optical floating zone method and cut to a disc 
shape with both surfaces polished parallel and to a mir¬ 
ror finish. The electric field profile of the THz pules after 
it transmits though the sample and an empty metal aper¬ 
ture are recorded. The ratio of the Fourier transforms of 
these traces can be further analyzed to give the complex 
transmission function of the sample, upon which further 
analyses are carried out. 

Microwave experiments are performed by the cavity 
perturbation technique, whereupon the resonance char¬ 
acteristics (center frequency and quality factor Q) of a 
superconducting resonator cavity are changed by the in¬ 
troduction of the sample. These changes are related to 
the complex generalized susceptibility of the sample un¬ 
der test. In both these techniques the complex suscepti¬ 
bility is measured, which gives both the reactive and the 
dissipative response of the material. Further details of 
these experiments can be found in the SI. 
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Supplementary Information 

A. Time-domain THz spectroscopy data analysis for both 1 and P 7^ 1 

In a time-domain THz spectroscopy (TDTS) measurement with sample in vacuum, the transmission (neglecting 
multiple reflections inside the sample) is given by the expression 


t = 


^Z 0 Z (n 

TNTW exp{ '~ 



(i) 


Here Z = is the wave impedance of the sample, Z 0 is the impedance of free space and n = ^f\xe is the complex 
index of refraction of the sample. For magnetic insulators, both the dielectric constant and magnetic susceptibility 
(Xm = yj) may show nontrivial features in the THz range. Therefore the commonly used assumption in the analysis 
of optical spectroscopy data that p, = 1 does not apply here. Disregarding this possibility could lead to misassignment 
of magnetic features as aspects in the charge conductivity ©■ Here we present a complete treatment of the situation. 

Expanding the above equation, we get 


t = 


_ / ,W e l x ^ l ) ujd \ 

> expyi ). 


(\Ie + a / m ) 2 


( 2 ) 


With n = 1 + Xrn = 1 + Xm + an d the reasonable assumption of Xm 1, this can be rewritten as 


+ VeU+Xm/ 2 ) (. (v^C 1 Xm/2) l)icd^ ^ 

(V~e+l + Xm/2) 2 P{ c > {S) 

Note that e and Xm are both complex quantities so that both may contribute to optical dissipation and phase 
retardation. From this equation it is obvious that the principle contribution of the optical response will typically 
come from the exponential dependence. 

As discussed in the main text, a separation of energy and temperature scales for the spin and lattice degrees of 
freedom allows us to extract the complex magnetic susceptibility from the TDTS measurements for materials like 
Yb2Ti20y. In the present case, the lowest infrared active optical phonon is located around 2.25 THz and so it is 
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FIG. 1 Effective index of refraction as a function of temperature Effective real index of refraction obtained from TDTS 
measurement, when setting p, = 1 in the data analysis. This data describes the total reactive response from the sample. Note 
the deviation from the temperature dependent trend below 15K. The deviation is more pronounced for lower frequencies. 


reasonable that the lattice becomes inert for temperatures much lower than this energy scale (~ 100 K) Therefore 
one can assume that the observed low temperature dissipation has a magnetic origin. This assumption is supported by 
temperature dependence of both the amplitude as well as the phase of the complex transmission function t measured 
in the TDTS experiments, as shown in Fig. 2(a) of the main text and the Supp. Fig. 1 here. Here we plot in Fig. 1 
the ‘effective’ index of refraction obtained from TDTS measurement by setting // = 1 in the data analysis. A close 
inspection of the data reveals that the temperature dependence of the reactive response changes sign at a temperature 
around 18 K, which is again indicative of the onset of magnetic degrees of freedom at this temperature range. 

With this information, then the temperature dependent changes to the optical response at temperatures lower than 
18 K can be attributed to a finite value of \ m in the THz frequency range, instead of a charge response. To isolate 
Xm from the complex transmission function, we choose a spectrum at T re f = 18 K as a reference spectrum. In the 
following data analysis, we assume that 18 K is the onset of magnetic response in the spectra, and the dielectric 
response does not change appreciably for lower temperatures. As discussed above, this assumption is supported by 
the energies of the optical phonons as well as the features of the TDTS data. 

Following the above discussion, we can extract y m (w) from the complex transmission via the expression 


ln( 


t(T ) 

t(T re f) 


) =i 


n'xmud 

2c 


(4) 


Here n! is the real part of the index of refraction of the sample obtained at the reference temperature. With this 
scheme of data analysis, the magnetic response can be obtained for quantum magnetic materials in suitable systems 
with similar separation of temperature and energy scales for dielectric and magnetic responses. Here we would like 
to emphasize that aside from this assumption the treatment presented here is completely model independent, thus it 
applies to all suitable materials and does not depend on the system being in a particular ground state. 


B. Microwave cavity data analysis 

Microwave cavity measurements were performed on several Yb 2 Ti 207 samples to measure the complex magnetic 
susceptibility at a frequency lower than that available in the THz. As discussed in the main text, these measurements 
were useful in constraining the Kramers-Kroning consistent fits to the complex susceptibility. In a typical experiment, 
one measures the shift in the complex resonant frequency Aw = Aw — i of a microwave cavity upon the introduction 
of a sample ©■ Here Aw = w s — wo is the shift in the cavity’s resonant frequency and AT = T s — To is the shift in 
the cavity’s bandwidth, or the full width at half maximum of the resonant peak. The subscripts “s” and “0” denote 
cavity with a sample and empty cavity respectively. If the introduction of the sample is a small perturbation on the 
cavity’s electromagnetic fields (e.g. the so-called “adiabatic” regime), it can be shown that the complex frequency 
shift is directly related to a complex generalized electromagnetic susceptibility ( of the sample, 

= -7 gC 


Aw 

w 0 


(5) 
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where j g is the geometrical factor that depends on the resonance mode of operation and the ratio of sample to cavity 
volume. Although the sample’s size is usually quite small for such measurements, a high quality factor of the cavity 
Q = ^ ss 10 4 - 10 s , allows for an exceptional level of sensitivity. For details of the method see Ref. [J, 

In our experiments, we use a superconducting NbTi cylindrical cavity, with T c ss 8.5 K, wo/27t « 18.5 GHz, Q 
~ 10 5 , designed to resonate in the TEon mode. The sample is mounted at the magnetic field anti-node at the center 
of the cavity where the electric field is zero, sitting atop a sapphire rod via a small amount of thermal grease. The 
cavity sits at the bottom of a He 3 cryostat capable of reaching a base temperature of 500 mK. The generation and 
detection of the microwave signals are performed by a network analyzer. The microwave response from the cavity is 
returned and fit to a Lorentzian function to extract the temperature dependent resonant frequency and bandwidth. 

In general, the susceptibility (£) of Eq. [5] is a generalized complex electromagnetic susceptibility which has contri¬ 
butions from both the sample’s electric (Xe) and magnetic (\ m ) susceptibilities the degree of which depends on the 
sample’s position within the cavity. In this regard it is important to note that, £ contains an electric susceptibility 
contribution even when the sample is placed in a pure magnetic field and vice versa. As with the TDTS experiments, 
both contributions must be taken into account. For small and insulating samples such as Yb 2 Ti 2 C> 7 , the electromag¬ 
netic fields within the cavity penetrate the sample’s volume completely and fc 0 a <C 1, where A : 0 = — is the microwave 
wavevector in vacuum and a is a characteristic sample dimension. It can be shown (|3]) that in this limit (known as 
the depolarization regime) and at a magnetic field maximum, the sample’s electromagnetic susceptibility can be well 
approximated by, 


c 


Xm 

1 Y 


1 Xm + 1 

10 (1 + N mX mf 


(k 0 a) 2 Xe 


( 6 ) 


where N m is the demagnetization factor of the sample that depends on sample shape. Thus in certain limits, electric 
and magnetic effects may be isolated from each other. 

For the Yb 2 Ti 2 07 sample investigated here, we estimate that k 0 a ~ 0.2. Similar to the TDTS analysis discussed 
above, we assume that changes to the electric and magnetic susceptibilities dominantly occur in two distinct temper¬ 
ature ranges, T > 18 K and T < 18 K respectively. However, our conclusions are not particularly sensitive to the 
choice of temperature scale, given that majority of the magnetic susceptibility signal onsets below 10 K. For T > 18 
K, assuming Xm = 0, plugging Eq. E into Eq. [5] gives, 


— = -Tx(fco a) 2 Xe, T > 18 K (7) 

UJq IU 

thus the frequency shift in this temperature range results from the sample’s (primarily real) electric susceptibility 
(polarizability). Because Yb 2 Ti 2 07 is a good electrical insulator, without any dielectric anomalies at low temperatures, 
we can assume that below T = 18 K the real part of Xe is constant while the imaginary part is zero. In the present 
case, because 1 and n ~ ^ the second term in Eq. [ojis approximately constant below 18 K, giving 


-= — 7 „--A?-[-constant, T < 18 K. ( 8 ) 

^0 1 H - NmX.m 

Thus, Xm can be extracted using Eq. [ 8 ]if Xe, 7 g, and N m are all known. In our case we extracted the electric suscep¬ 
tibility by linearly extrapolating our THz data into the low frequency regime (x e of 52 and 66 at room temperature 
and 18 K respectively). There is about a 17% rise in x' e from the microwave range to 1 THz. We then used the shift 
in resonant frequency from room temperature to 18K of the loaded cavity and relate a change in Aw to a change in 
Xe over this range to find 7 g , 


A(—) = -^(k 0 a) 2 A Xe . (9) 

wo 10 

The absolute shift between the unloaded and loaded cavity was not used in order to reduce the error caused by a 
resonant frequency offset that results from removal and replacement of the cavity’s bottom plate when inserting a 
sample. The demagnetization factors of our samples were estimated from the sample’s dimensions. An extrapolation 
in frequency of Xe from the THz range to the microwave range should be valid in a good electrical insulator like 
Yb 2 Ti 2 C> 7 . As it was assumed that the magnetic susceptibility of the sample was negligible above 18 K, the magnetic 
susceptibility below 18 K could then be calculated from Eq. [ 8 ] when using the relative shift between the 18 K resonant 
frequency and the low temperature value. 
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Temperature (K) Temperature (K) 

FIG. 2 Temperature dependence of the microwave cavity data Temperature dependence of the real (a) and imaginary part (b) 
of the complex frequency shift for sample A (red) and sample B (blue), as well as the temperature dependence of the extracted 
real (c) and imaginary (d) magnetic susceptibility. 


Sample 

Dimensions (mm) 

Geo. Factor y 9 

Demag. Factor N m 

A 

0.95 x 1.015 x 0.582 

5.6xl0" 4 

0.49 

B 

0.95 x 1.015 x 0.479 

4.9x 10 -4 

0.52 


TABLE I Dimensions, extracted geometrical factors, and calculated demagnetization factors of the two Yb 2 Ti 2 C >7 samples 
measured (see text for details.) 


To ensure accuracy of the microwave measurements, two Yb 2 Ti 2 C >7 samples with different thickness were measured 
and the results compared to each other. Table |T] summarizes the dimensions, measured geometrical factors, and 
calculated demagnetization factors of the two samples. After measuring, sample A was polished to a thickness of 
0.479 mm and remeasured as sample B. Figs. SI2 (a) and (b) show the temperature dependence of the change of 
the resonant frequency and bandwidth from the two samples respectively. The superconducting transition of the 
cavity is seen at T « 8.5 K, above which there is significantly reduced resolution. The reduced resolution above the 
cavity’s transition does not significantly effect our extracted low temperature susceptibility as the signal was small 
above 8.5 K. Strong signatures from the sample’s magnetic correlations are seen below T sa 8 K in both the resonant 
frequency and bandwidth of the loaded cavity. The difference in magnitudes of the shifts from the two samples are 
caused by different geometrical and demagnetization factors. However, the general shape of the curves and even 
more importantly, the direction of the shift is consistent between measurements. These data and the parameters 
summarized in table [I] were then used to extract the magnetic susceptibility of each sample. 

In Fig. SI2(c) and (d) we plot the extracted real and imaginary magnetic susceptibility. The data shows a general 
trend of increase upon cooling, which is what one expects for this material. We found less than a 10% difference in 
the absolute magnetic susceptibilities between the two samples in the temperature range 1.5 K < T < 8 K, which 
should be considered quite good. The difference could stem from uncertainties in the extracted geometrical factors or 
in the demagnetization factors resulting from approximating our rectangular prism samples as ellipsoids. 
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Frequency (THz) 

FIG. 3 Real part of \m measured form magneto-optical TDTS Real part of obtained from a TDTS measurement from 
Yb 2 Ti 207 in field, with temperature at 1.6 K with H ac || H dc -L ^thz (Voigt geometry). With the application of an external 
field, the zero crossing in x'm is pushed to higher frequencies. 


C. Magneto-optical TDTS 

To further investigate the reactive response of the monopole conductivity, we performed a magneto-optical TDTS 
measurement with H ac |j H ^ _L k thz (Voigt geometry). The measurement was conducted with a superconducting 
magnet with optical access. The reactive and dissipative response from the measurement at 1.6 K are shown in Fig. 
SI3 and Fig. 2(b) of the main text, respectively. Similar data can be found in Ref. ({5]). 

With the application of an external magnetic field, the diffuse zero field magnetic spectral weight in the dissipative 
response becomes sharper and is pushed to higher frequencies. Corresponding features are also seen in the reactive 
response. Here similar data analysis as discussed in the first part of the supplementary information is applied to 
extract Xm in the magneto-optical measurements. As clearly shown in the figure, with the application of magnetic 
field, the spectra features smoothly shift to higher frequencies with the zero crossing in x'm move to higher frequencies 
where reliable THz data are available. The data provides strong evidence for the magnetic origin of the main features 
as well as showing a sign change in the reactive magnetic response. 


D. Extended Ryzhkin model 

In analogy with previous work done on proton disorder and dielectric relaxation in water ice, Ryzhkin (|B]) derived 
an expression for the magnetic relaxation through monopole motion in classical spin ice. Starting from the formalism 
for the local entropy production in irreversible process he derived the monopole flux to be 


J z — M^z,m(9z,mH M (10) 

where i refers to positive and negative monopole charge species, /i is the monopole mobility, n.; >m is the monopole 
density, and rj z = ±1. q i m is the monopole magnetic charge which is set by the condition that the spin dipole moment 
fi m equals q m a/2 where a is spin-spin bond length. Ryzhkin’s expression differs from a simple transport model of free 
charge by the inclusion of the second term which includes a dependence on a configuration vector f2 that is related to 
the system magnetization as fl = M /q m . Finite magnetization reduces entropy and provides a thermodynamic force 
that opposes the current. This reaction force originates in the configurational entropy of the monopole vacuum and 
prohibits a true dc current even in the absence of sample boundaries. *!> is a constant of proportionality derived to be 
^akgT in the context of water ice models that retains its relevance here 0- As magnetic relaxation proceeds in a 
spin-ice through monopole motion the configuration vector is related to the history of the monopole current as 






12 


fl(t) - fi( 0 ) 



( 11 ) 


Substituting the Fourier transform of Eq. ED into Eq. [TO] yields a Debye-like relaxation form for the susceptibility 
{Xm = X'm + iXm) 


Xm(w) 


1 — ioJT 


( 12 ) 


with a relaxation time 




in which n m is the total (thermally excited) monopole density n+ m + n_ m . As 
discussed in the main text, the leading dependence of the susceptibility in powers of 1 /(—iu) is a diagnostic for 
the dominant term in the equations of motion. Accordingly Eq. 12 is characterized by a 1 /(—iui) fall-off at high 
frequencies that reflects with a dominant dissipative response. At low frequencies \ m is constant and real showing 
(e.g. the zeroth power of 1 /(—iui)) that the dominant effect in the dc limit comes from the reaction force. Note that 
both real and imaginary parts of Eq. 12 will be positive for all frequencies e.g. there is no sign change in x'm- 

Although such Debye-like relaxation has been found to be a good description of the relaxation processes present in 
classical spin ice (Ho 2 Ti 2 07 and Dy 2 Ti 2 C >7 ) at the low (kHz) frequencies used in most experiments, it is not strictly 
speaking a fully mathematically consistent response function as it does not fall off fast enough at high ui to satisfy, 
for instance, the zero temperature first moment sum rule f 0 °° dui u>Xmaa( q- w ) = [H, 5“ q ]]) for local spins 

(SiU). Here 5“ = ^ ■ e ® q ' Rj '5“ is the Fourier transform of lattice spin operators and T~L is a spin hamiltonian. The 
evaluation of the first moment sum rule depends on the details of the exchange interaction, but being equivalent to a 
double commutator of the spin hamiltonian with spin operators, will yield some finite value. 

A similar violation of the sum rule is a well-known defect of the Debye model in its application to electronic charge 
relaxation in dielectrics m- In general, the first moment of the electric susceptibility must obey the sum rule 


J 0 °° dui tux" = § where n e is the total charge density and m e the electron mass. Similar to the Ryzhkin model, 
the Debye model does not fall off fast enough at high u> to satisfy the sum rule. The addition of inertial terms to the 
Debye model in the classical equations of motion does give a Xm a high frequency asymptote that goes like l/(— iui ) 2 
which insures its first moment’s integrability (jlO|) . Although at low frequencies in strongly dissipative media, inertia 
can be neglected, at high enough frequencies it must become relevant to satisfy the sum rule. Frequently in optical 
studies of charge systems a restricted low energy sum rule is applicable, where the upper limit of integration is taken 
to be finite, but then n e is replaced by the density of electrons in the low energy sector, and m e is replaced by a mass 
that is renormalized by band and/or interaction effects. 

Therefore, for spin ice, irrespective of their exact form, effects beyond Ryzhkin’s treatment must become relevant 
at high enough frequency to satisfy the first moment sum rule. We have argued in the main text that inertial effects 
must be included to understand magnetic relaxation in the quantum spin ices. In a similar fashion to charge in 
dielectrics, the Ryzhkin expression Eq. [12] can be amended by including a phenomenological inertial term to Eq. |10| 
As discussed in the main text, it will have the form of the final term in the expression 


J i — Jz/7 


(13) 


where 7 is a current relaxation rate. 

From the standard definitions, q m /H = 7 m. In what follows, m will be an effective inertial mass which arises in 
the low energy sector through non-Ising exchanges that lead to monopole tunneling. There is no more contradiction 
inherent in using a classical model to describe a quantum spin ice than there is to use a classical model to describe 
the conduction of thermally excited charge in a semiconductor where classical inertia also arises through inherently 
quantum tunneling. Solving in the same fashion as in the above Ryzhkin case gives a classical equation of motion, 
the terms of which are instantly familiar and map to the form of a damped harmonic oscillator. 

M + 7 M+—M = ^H. (14) 

m m 

Solving for the susceptibility we have 


_ dm/® _ 

1 — iio'ym /n m *I> — u] 2 m/n m & 


XmM 


(15) 
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which in terms of r (and 7 ) reads 


Xm(w) = 


q 2 J* 


1 — iojT — lc 2 t jx 


With appropriate substitutions (wg = 


rim& 


and Wp = 


Eq. [l5] 


(16) 


can be seen to be equivalent to the Drude- 


Lorentz equations that describe the response of a classical electric charge oscillator. 


XmM 


LOt, 


Wg — ui 2 — iui^f 


(17) 


In the limit where 7 -7 00 Eq. [16] recovers Ryzhkin’s expression for Debye-like relaxation. In this limit the zero 
crossing in x" will move to infinity. 

From Eq. [16] for the susceptibility and in analogy with charge conductivity we can define a magnetic monopole 
conductivity k(ui) = —iuiXm((o) = k' + in". As discussed in the main text, the quantity coXm. is the quantity directly 
measured in TDTS experiments. Due to mapping of our expressions Eqs. |15| and [l6| to the Drude-Lorentz model, the 
magnetic conductivity must obey a low energy sum rule where 


r 7 r 

/ du> K {bS) = — 
Jo 2 


7T 


m 


(18) 


Here ui* is a energy cut-off that must be smaller than a scale on the order of the Ising exchange parameter J zz , but 
much larger than 7 . Note that this is extension to the usual first moment sum rule for spin systems as k = —icoxm- It 
is also important to note that this monopole conductivity k is not the same as the spin conductivity which is defined 
as the response of a spin current to the magnetic field gradient in the small momentum limit mm- 

In Fig. SI4 we plot the results of function Eq. 16 for both susceptibility Xm and the magnetic conductivity k for the 
range of r and 7 we believe is relevant to Yb 2 Ti 20 y. Here we use example parameters 1/r = 4 and 7 = 16 on both 
linear, log, and log-log scales. These plots highlight the utility of analyzing both of these quantities. We see that the 
dissipative part of the susceptibility x" peaks near 1/r, which is the usual Debye-like relaxation behavior. In Xm the 
scale of 7 appears only as a subtle change in the high frequency power law, which manifests itself as a change of slope 
on the log-log plot. An experimentalist doing typical low frequency susceptibility or neutron scattering experiments 
(that measure only would be likely to be completely unaware of the high frequency scale where inertial effects 
become relevant. The existence of inertial effects are manifest due to the sign change in x' m (which are absent in 
the Debye-like Ryzhkin expression). The magnetic conductivity k exhibits both frequency scales prominently. The 
imaginary part of n exhibits a negative extremum at 1 /r. 7 is seen as the frequency where k' > k" and n" exhibits a 
maximum. The upper right panel of n on linear scale demonstrates that if one is measuring at frequencies well above 
1 It the response will be indistinguishable from a Drude-like transport and appear as a monopole metal. In that case 
the width of the peak in the real part of the monopole conductivity would give 7 . 

The lower left panel also demonstrates the point made in the text that the bounding behavior of Xm is a power 
law in 1 /(i(o), where the power indicates what term in the classical equation of motion is dominant at a particular 
frequency. This is because there are three “response” terms in the equation of motion, which leads to three terms in 
the denominator of the susceptibility t hat can be written as increasing powers in iui. In different frequency regimes, 

dominate. Distinct frequency regions of l/(iw) 0 , 


different terms on the bottom of Eq. 


16 


l/(iw), and 1 /(ilo) 2 are 


therefore apparent in the plots of Fig. SI4. 

A few more comments are in order about the applicability of our model for % m . Our approach above was to first 
define macroscopic variables (J, M etc.) and then combine them in a way that was consistent with both an entropic 
restoring force and inertial effects. The perceptive reader may make the reasonable objection to this approach that 
the long time averaging which is necessary to define an entropic force is inconsistent with the short time scales on 
which inertial effects will become apparent. An alternative and more rigorous approach would be to first write down 
a Langevin type equation that includes inertia and then do the averaging in the Langevin style. In contrast our 
approach was one where the order of these operations was essentially reversed. It is not a priori clear if the steps of 
including the dynamics as such and then averaging or vice versa commute. In general they do not. 

Of course, this issue is not unique to the present case. Such considerations are ubiquitous when a high frequency 
response is considered in the presence of inertial forces (e.g. polymer relaxation, high frequency phonon propagation, 
etc.). For instance, related issues arise in the inclusion of inertial effects in the context of dielectric relaxation. A 
rigorous Langevin treatment of an assembly of non-interacting fixed axis rotators with finite moments of inertia yields 
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FIG. 4 Simulations of frequency dependent monopole susceptibility and conductivity Plots of the susceptibility \m and the 

magnetic conductivity k with example parameters 1/r = 4 and 7 = 16 on both linear, log, and log-log scales. 


a susceptibility that has the complicated form of a infinite set of continued fractions (13) (Similar expressions were 
obtained earlier by different methods by Gross m and by Sack (1 161) 1. However, for small inertia, this expression 
can be approximated by the first convergent, and is known as the Rocard equation, after Rocard who in 1933 first 
derived it da when including inertial effects in the context of the Debye theory. The Rocard equation is completely 
equivalent to the expression we have used above. This example shows that there are cases and limits in which the 
ordering in the treatment of including inertia and ensemble averaging does not effect the end result. Unfortunately, 
it is not known at this time, to what extent these effects matter or how to perform a Langevin style analysis for 
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FIG. 5 Temperature dependence of monopole mass and mobility Temperature dependence of (a) the effective monopole mass 
in units of electron mass, and (b) the monopole mobility. 


dynamics in the spin ices. This is a very open area for theoretical inquiry. 


E. Monopole mass and mobility 

Following the discussion in the last section, the effective mass and mobility of magnetic monopoles can be obtained 
from the fitting parameters of the monopole conductivity. The monopole mass can be estimated from the values of 
the magnetic spectral weight, the monopole charge, as well as the numerical results for the monopole density. The 
result, in unit of the free electron mass, is plotted in Fig. SI5(a). As discussed in the main text, at low temperatures 
the values shown here are in good agreement with a low temperature monopole “band mass” of approximately 2200 
m e calculated in a model with charges hopping on diamond lattice with amplitude J z ± and using the parameters of 
Ref. (|T7|) . However, at higher temperatures the effective mass increases strongly and becomes of order 10,000 m e 
at 10K. We believe that these large values do not have any physical significance and is an indication that one is no 
longer in a temperature regime where a monopole model has validity. At temperature above the mean field T c the 
fitting parameters should be considered phenomenological only. 

The monopole mobility can be obtained from /i = q m /jm once the effective mass is available. The results are 
shown in Fig. SI5(b). A value several orders of magnitude larger than the ones found in classical spin ice is obtained 
here m- The huge contrast of monopole mobilities in two different classes of spin ice materials point to essential 
different physics that governs the low temperature dynamics of the monopole excitations. As discussed in the main 
text, although we have performed these fits up to temperature of order 10K, we believe these quantities should only 
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be taken to describe actual physical monopole parameters below the nominal mean-field transition temperature scale 
of ~ 4K. 
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